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Abstract 



An 50(10) grand unified model considered previously by the authors fea- 
turing lopsided down quark and charged lepton mass matrices is successfully 
predictive and requires that the lightest two right-handed Major ana neutri- 
nons be nearly degenerate in order to obtain the LMA solar neutrino solu- 
tion. Here we use this model to test its predictions for baryogenesis through 
resonant-enhanced leptogenesis. With the conventional type I seesaw mecha- 
nism, the best predictions for baryogenesis appear to fall a factor of three short 
of the observed value. However, with a proposed type III seesaw mechanism 
leading to three pairs of massive pseudo-Dirac neutrinos, resonant leptogen- 
esis is decoupled from the neutrino mass and mixing issues with successful 
baryogenesis easily obtained. 
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I. INTRODUCTION 



The accumulation of refined data on both neutrino oscillation and quark flavor physics 
has posed increasingly severe tests for quark and lepton models which attempt to explain the 
mass and mixing data. This is especially true for Grand Unified models which must relate 
both the quark and lepton sectors. In fact, at present several such models [1] do survive these 
tests given the present levels of experimental precision for the Cabibbo-Kobayashi-Maskawa 
(CKM) quark and Maki-Nakagawa-Sakata (MNS) lepton mixing matrix elements. Two 
other hurdles for these Grand Unified models involve the issues of flavor-changing neutral 
currents, for example in \i — > e + 7 and r — > /x + 7 decays, and the survival of baryogenesis 
in the early Universe. It is the latter hurdle which we wish to address in this paper. 

It appeared early on that Grand Unified Theories (GUTs) provided the necessary condi- 
tions for successful baryogenesis as postulated by Sakharov [2]: baryon- violating interactions 
which involve both C- and CP- violation and which occur out-of-equilibrium in the early Uni- 
verse. Since then we have learned that the baryon-violating interactions which occur near 
the GUT scale conserve B — L, so that any net baryon number generated by them can 
be washed out by sphaleron B + L interactions which occur in thermal equilibrium with 
the expanding universe. Only if the net baryon number is generated with A(B — L) 7^ 
interactions will it not be erased by the sphaleron interactions [3]. 

For this reason, Fukugita and Yanagida suggested that leptogenesis may play a necessary 
primary role for baryogenesis [4]. An excess in lepton number generated by the lepton- 
violating decays of the heavy right-handed neutrinos can be converted into a baryon excess 
by sphaleron interactions in thermal equilibrium above the critical electroweak symmetry- 
breaking temperature. In this decay or direct "e"' CP-violating scenario, the CP violation 
is generated through an interference between the decay tree graph and the absorptive part of 
the one-loop decay vertex. This mechanism requires very heavy iVj, % — 1, 2, 3 right-handed 
neutrinos, with the lightest being heavier than roughly 10 9 GeV. Thus with a moderate 
hierarchy of heavy right-handed neutrinos and with the lightest satisfying the above bound, 
satisfactory leptogenesis can be generated. 

However, it has been observed that successful leptogenesis may also arise through indirect 
"e" CP-violating mass-mixing effects in the decays of two quasi-degenerate heavy right- 
handed neutrinos [5]. In other words, resonant enhancement of the lepton asymmetry can 
be generated if the two neutrinos are sufficiently close in mass that level crossing can occur. 
Ellis, Raidal, and Yanagida as well as Akhmedov, Frigerio and Smirnov have performed 
phenomenological studies [6] of neutrino mass matrices and shown that sufficient lepton 
asymmetry can be generated in the case that the level crossing involves the two lighter 
right-handed neutrinos, N\ and N 2 ; moreover, the quasi-degenerate masses can lie as low as 
10 8 GeV, at least a factor of 10 below the bound obtained in the generic hierarchical case. 
Recently Pilaftsis and Underwood have proposed a model of the lepton sector in which the 
nearly degenerate neutrino pair is as light as 1 TeV [7]. 

It is important to test this suggested scenario in a realistic model of quark and lepton 
masses in as quantitative a manner as possible. At issue is whether or not the desired results 
for both leptogenesis and neutrino oscillations can be obtained simultaneously for the model 
in question. In this paper we explore leptogenesis with such resonant enhancement in an 
SO (10) GUT model proposed several years ago by us in collaboration with Babu [8] with 



2 



a series of refinements by us since then [9]. This model is numerically very predictive and 
is found to explain accurately the present data on quark masses and mixings as well as 
the large mixing angle (LMA) Mikheyev-Smirnov-Wolfenstein (MSW) [10] solar neutrino 
solution with near maximal atmospheric neutrino mixing. The lopsided texture of the 2-3 
submatrices of the down quark and charged lepton mass matrices neatly explains the small 
mixing in that sector of the CKM matrix, i.e., small V c b and Vt. s , as well as the near maximal 
mixing in that sector of the MNS matrix [11], i.e., large U^. The LMA solar neutrino solution 
arises from the type I seesaw mechanism as an interplay between the Dirac neutrino matrix 
and the right-handed Majorana mass matrix, M R . In fact, the 2-3 submatrix of M R has 
a zero determinant which nearly replicates the corresponding sector for the right-handed 
Dirac neutrino matrix product, M^M^. With the full structure of Mr spelled out to give 
the LMA solution, one finds the lighter two masses N\ and N 2 are nearly degenerate, so at 
least the possibility of resonant enhancement of leptogenesis exists in this model. We shall 
examine how successful the model can be in achieving leptogenesis sufficient to generate the 
desired amount of baryogenesis observed in Nature. 

As an alternative approach, we can reformulate the model with the type III seesaw 
mechanism developed in [12,13]. Whereas in the usual type I see-saw mechanism there 
are three superheavy Majorana neutrinos, in the type III see-saw, there are six superheavy 
neutrinos that typically form three pseudo-Dirac pairs. The desired resonant enhancement 
of leptogenesis can result from the mixing of the two neutrinos of the lightest pseudo-Dirac 
pair. In this scheme, the LMA solar neutrino mixing solution emerges rather naturally 
without fine-tuning of the Dirac and Majorana neutrino matrices, and the heavy lepton 
decay asymmetry is effectively decoupled from the light neutrino neutrino mass and mixing 
issues. As a result we shall see that satisfactory leptogenesis is easily obtained in this 
scenario. 

In Sect. II we summarize briefly the features of the SO (10) model and take the oppor- 
tunity to update the input parameters so as to give even better agreement with the quark 
and lepton mass and mixing data. The formulas relevant for leptogenesis involving two 
quasi-degenerate right-handed neutrinos are summarized in Sect. III. We then apply these 
formulas to the model, modify the model slightly and compare the leptogenesis and oscilla- 
tion results with experimental information. In Sect. IV the necessary formulas are presented 
for the type III seesaw and numerically satisfactory examples are given. Conclusions follow 
in Sect. V. 

II. 50(10) MODEL WITH 17(1) x Z 2 x Z 2 FLAVOR SYMMETRY 

The GUT model in question [8,9] is based on the grand unified group 5*0(10) with 
a U(l) x Z 2 x Z 2 flavor symmetry. The model involves a minimum set of Higgs fields 
which solves the doublet-triplet splitting problem. This requires just one 45 h whose VEV 
points in the B — L direction, two pairs of 16//, 16j/'s which stabilize the solution, along 
with several Higgs fields in the 10# representations and Higgs singlets [14]. The Higgs 
superpotential exhibits the U(l) x Z 2 x Z 2 symmetry which is used for the flavor symmetry 
of the GUT model. The combination of VEVs, (45 h )b-l, (1(16 h )> and (1(16 H )) break 
50(10) to the Standard Model. The electroweak VEVs arise from the combinations v u = 
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(5(10h)) and Vd = (5(10h)) cos 7 + (5(16' H )) sin 7, while the combination orthogonal to Vd 
gets massive at the GUT scale. As such, Yukawa coupling unification can be achieved at 
the GUT scale with tan/3 ~ 2 — 55, depending upon the 5(10h) — 5(16h) mixing present 
for the Vd VEV. In addition, matter superfields appear in the following representations: 
I61, 16 2 , 16 3 ; 16, 16, 16', 16 7 , 10i, 10 2 , and l's, where all but the 16 ; (i = 1,2,3) get 
superheavy. 

The mass matrices then follow from Froggatt-Nielsen diagrams [15] in which the super- 
heavy fields 16, 16, 16', 16', 10i, 10 2 , and 1 are integrated out. These diagrams, cf. [16], 
correspond to the following effective Yukawa operators for the indicated Dirac mass matrix 
elements: 
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where the subscripts to the brackets indicate how the fields inside the brackets are contracted. 
The Dirac mass matrices for the up quarks, down quarks, neutrinos and charged leptons are 
then found to be 
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The above textures give the Georgi-Jarlskog relations [17] between the quark and lepton 



GUT scale masses, 



ml/3, 



3m° with Yukawa coupling unification holding for 



tan/3 ~ 5. (Here the Dirac matrices are written with the convention that the left-handed 
fields label the rows and the left-handed conjugate fields label the columns. The opposite 
convention was used in some earlier references for this model. Hence the matrices here are 
the transpose of those given in those earlier papers.) 

All nine quark and charged lepton masses, plus the three CKM angles and CP phase, 
are well-fitted with the eight input parameters (the older choice is indicated in parentheses) 



mu ~ 113 GeV, 

a = 1.83 (1.78), 

5 = 0.00946 (0.0086), 

= 119.4° (126°), 



m D ~ 1 GeV, 

e = 0.147 (0.145), 
5' = 0.00827 (0.0079), 
r;=6x 10~ 6 (8 x 10- 6 ), 



(3) 



defined at the GUT scale to fit the low scale observables after evolution downward from 
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A better overall agreement with experiment [18] is obtained with the new parameter values 
aside from the electron mass which is most sensitive to small corrections. With no extra 
phases present other than the one appearing in the CKM mixing matrix, the vertex of the 
CKM unitary triangle occurs near the center of the presently allowed region with sin 2/3 ~ 
0.72. The Hermitian matrices M v M\j, M d M ] d , and M N M ] N are diagonalized with small 
left-handed rotations, while MlM\ is diagonalized by a large left-handed rotation. This 
accounts for the small value of V c b = (Uu L U DL ) cb , while \U^\ = KU^U^)^ will turn out 
to be large for any reasonable right-handed Majorana mass matrix, Mr [8]. 

The effective light neutrino mass matrix, M v , is obtained from the type I seesaw mecha- 
nism once the right-handed Majorana mass matrix, Mr, is specified. The large atmospheric 
neutrino mixing <-> v T arises primarily from the structure of the charged lepton mass 
matrix Ml, while the structure of the right-handed Majorana mass matrix Mr determines 
the type of u e v^, v T solar neutrino mixing, so that the solar and atmospheric mixings 
are essentially decoupled in the model. The LMA solar neutrino solution is obtained with 
a special form of Mr, as will be seen in a moment. However, this special form can be 
explained by the structure of the Froggatt-Nielsen diagrams [16]. The most general form 
for the right-handed Majorana mass matrix considered in [16] which gives the large mixing 
angle (LMA) solar neutrino solution is 

(c 2 rf —her] arj \ 
-berj e 2 -e A R , (5) 
arj —e 1 J 

where the parameters e and rj are those introduced in Eq. (2) for the Dirac sector. With 
a 7^ b = c, for example, the structure of Mr arises in the following way. The VEV of one 
particular Higgs singlet that has AL = 2 contributes to all nine matrix elements giving a 
factorized rank 1 form. The VEV of a second Higgs singlet also breaks lepton number, but 
contributes only to the 13 and 31 elements of Mr. 

Given the right-handed Majorana mass matrix above, the seesaw formula results in 

M y = -Mnmj-mj; 



As a numerical example, with just three additional input parameters: a — 1, b = c = 2 and 
A R = 2.65 x 10 14 GeV, one obtains 
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leading to 



5.1 x 10" 3 eV, m 2 = 9.1 x 10 
M 2 ~ 2.3 x 10 8 GeV, 
= 2.6 x 10~ 3 eV 2 , 
_ = 5.6 x lO" 5 eV 2 , 
U e3 = -0.0172 - 0.00094i, 
J = 2.0 x 10~ 4 , 5 CP = 177°, 



m! = 
Mi ~ 
Am 2 2 

A»21 



-3 eV, m 3 = 52 x 10" 
M 3 = 2.7 x 10 14 GeV, 
sin 2 26 atm = 0.991, 
tan 2 0i2 = 0.48, 
sin 2 20 13 = 0.0012, 
Xi = -180°, X 2 = 90°, 



eV, 



(7) 



(8) 



to be compared with the present atmospheric, solar and reactor data and best-fit point in 
the LMA region [19] 



Am 2 2 ~ 2.6 x 10~ 3 eV 2 , sin 2 29 atm = 1.0 (> 0.92 at 90% c.L), 
Am 2 x = 7.1 x 10^ 5 eV 2 , tan 2 9 12 = 0.44. 



(9) 



In fact, the whole presently-allowed LMA region can be covered with a thin strip in the a — b 
plane given by b — c ~ 1.9 + 1.4a. In this region, sin 2 2#i 3 = 0.0004 — 0.0012. Although the 
prediction for U& is thus 3-5 times larger than the CKM mixing element V u b, a Neutrino 
Factory would be required to reach that range of sin 2 29 13 . 

It is interesting to remark that the hierarchy exhibited by the light left-handed neutrinos 
is a weak and normal one as is typical in 50(10) models. This is somewhat surprising, for 
both the Dirac and Majorana neutrino matrices, N and Mr, exhibit strong hierarchies. As 
is apparent these hierarchies nearly cancel each other in the type I seesaw mechanism. 

It is a simple matter to compute the effective mass parameter for neutrino-less double 
beta decay. For this purpose, we first note that the MNS neutrino mixing matrix is given 
by the product of the two unitary matrices which diagonalize the charged lepton and light 
left-handed neutrino mass matrices, i.e., Vmns = U\ j JJ Ul . The diagonalization occurs as 
follows: 

U T LL M L MlUl L = diag(m 2 , m 2 , m 2 ), 
U^ l M Vl U Vl = diag(mi, m 2 , m 3 ), 

where the mass eigenvalues are taken to be positive. Clearly an arbitrary phase trans- 
formation can be made on Ul l but not on U VL . Hence the mixing matrix Vmns has the 
form 

Vmns = U M ns®, $ = diag(exp % , exp % , 1), (11) 
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where Umns has the standard form with the Dirac CP phase appearing in the 13 element, 
while X\i X2 are the two Majorana CP phases. In terms of the above, the effective mass in 
neutrino-less double beta decay can then be written as 

(m ee ) = | Ej rrij {U M NS^)\j I 

= | J2j mjVj\UMNS,ij\ 2 \, 

where r]j = exp 2 * x -? for j = 1, 2 is the CP-parity of the jth lepton. For the example illustrated, 
Xi = —180°, X2 = 90° and (m ee ) = 0.57 meV, well below the present or future limit of 
observability, typical for a normal hierarchy spectrum. The lightest two neutrino states 
have opposite CP-parity. 

Note that the two lightest heavy right-handed neutrinos are nearly degenerate with 

(M 2 - Mi)/M 2 = 1.21 x 1(T 4 . (13) 

Actually they are nearly degenerate in magnitude only, their opposite relative signs that 
appear in the eigenvalues of Mr signify they have exactly opposite CP-parity, a feature which 
enables them to evolve downward from the GUT scale without their separation receiving 
large radiative corrections. The small relative separation of the quasi-degenerate states 
Ni and N 2 suggests that some resonance enhancement of leptogenesis may result from the 
level crossing. As we shall see in the next Section, this separation is too large by more 
than an order of magnitude to produce enough leptogenesis. Moreover, at least one of 
the four parameters a, b, c and Ar in the right-handed Majorana matrix defined in (5) 
must be complex in order to generate a lepton asymmetry. We shall attempt to modify 
the parameter assignments in M R to achieve satisfactory leptogenesis after discussing the 
required conditions in the next Section. 



III. LEPTOGENESIS WITH TYPE I SEESAW 

Here we present the basic formulas for calculation of the lepton asymmetry which results 
in baryogenesis with the type I seesaw mechanism. For this purpose, we can use as a guide 
the recent phenomenological studies carried out by the authors of ref. [6]. There is an 
important difference in our approaches, however. They chose to work in the charged lepton 
flavor basis with the right-handed Majorana mass matrix Mr diagonal, while our model 
was naturally developed in the 5*0(10) flavor basis. Hence we must transform our neutrino 
matrices and Mr to the basis in which Mr is diagonal in order to apply the conventional 
formulas in refs. [4-6]. 

The basic assumption is that a lepton asymmetry is generated by the CP-violating 
out-of-equilibrium decays of the right-handed neutrino iVj. Since the right-handed Majorana 
mass term violates lepton number by two units, Ni is identical to its charge conjugate state 
and can decay in two ways: 

Ni -> Vj + 0, Vj + f ; i,j = 1, 2, 3 (14) 

where is the Higgs field coupling the right-handed Ni to the light left-handed neutrino v r 
The Yukawa coupling involved is given by the ijth element of the Dirac matrix M' N in the 



7 



basis where the right-handed Majorana neutrino mass matrix M' R is diagonal. Let us label 
that Yukawa coupling by h'^. The lepton asymmetry is then traced to the imaginary part of 
the interference arising from the direct decay diagram and the one-loop diagrams depicted 
in Fig. 1: 
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When summed over both neutrino and charged lepton final states for all three families, the 
total lepton asymmetry is still given by the right-hand side of Eq. (15); hence we have 
omitted the subscript j on e*. The first term in f(x) arises from the self-energy correction 
while the second and third terms arise from the one-loop vertex. In the limit of quasi- 
degenerate neutrinos, f(x) for x = |M fc | 2 /|Mj| 2 reduces to 



lim/(a?) 









2(\M k 




\Mi\) 



(18) 



A resonance enhancement thus arises from the level crossing or near vanishing of the denom- 
inator of f(x) in the case of nearly degenerate masses. The lepton asymmetry e\ associated 
with the decay of Ni thus simplifies to 
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Mi 




16tt 


M 2 




Mi| 



X, 



21- 



(19) 



The lepton asymmetry e 2 associated with the decay of N 2 is found to be equal to that of e\ 
for nearly degenerate states. In either case, the enhancement is limited by the decay widths 
of the two states for ]|M 2 | - |Mi|| ~ I\/2 ~ T 2 /2, cf. [20], where 



— {^h'UMl 

07T 



(20) 



In Eq. (15) above, h! is the Yukawa coupling matrix for the Dirac neutrinos in the 
basis where the right-handed Majorana matrix M R is diagonal, which we shall call the 
primed basis. Denoting mass matrices written in that basis with primes, one has then 
h' = M' N I '{v sin /3) with v = 174 GeV. In Eq. (2) the mass matrices are given in the original 
50(10) flavor basis, which we call the unprimed basis. 

If the various neutrino mass matrices are diagonalized as follows 
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then one has that 



UJ, l M n U* Nr = diag(m u , 9m c , m t ), 

U^ L M u Uu L = diag(mi, m 2 , m 3 ), (21) 
Ul R M R U* MR = diag(M 1 , M 2 , M 3 ), 



M%M' N = U T M MiM N W M . (22) 



The transformation matrix Um r is uniquely determined provided we require the mass eigen- 
values in Eq. (21) be real and positive. The factor of 9 in the diagonalized Mfq can be 
understood by comparing and Mjj in Eq. (2). 

The lepton asymmetry produced by the decay of iVj is partially diluted by the lepton 
number- violating processes themselves [21]. This washout is determined as a function of an 
effective mass given by 

•fh i = {M%M' N ) ii l\M i \. (23) 
The washout factor /Cj for i — 1,2 with rhi in the range 1CT 2 — 10 3 eV is approximated by 

The lepton number asymmetry produced per unit entropy at temperature T > Mi ;2 , 
taking into account decays of both Ni and N 2 and their nearly equal washout factors, is 
then given by [22] 

n L 2k 1 € 1 g N T 3 

7T ( 25 ) 



S S IT 

90 g N 



2vr 4 g* 



Kiei. (26) 



We have used the expression for the entropy of the co-moving volume, s = (2/45)g*7r 2 T 3 . 
Here g^ = 2 refers to the two spin degrees of freedom of each decaying Majorana neutrino, 
while g* = 106.75 refers to the effective number of relativistic degrees of freedom contributing 
to the entropy in the absence of supersymmetric particles. Hence we find 

(I^)SM ^ 8M x 10 -3 Kiei _ (27) 

s 

The corresponding B — L asymmetry per unit entropy is just the negative of til/s, since 
baryon number is conserved in the right-handed Majorana neutrino decays. While B — L 
is conserved by the electroweak interactions following those decays, the sphaleron processes 
violate B + L conservation and convert the B — L asymmetry into a baryon asymmetry. 
Following the work of Harvey and Turner [3] the connection is 

ns „ 24 + 4N H n L 

s 66 + 13N H s ' 1 ; 
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where Nh is the number of Higgs doublets. Again in the absence of supersymmetric con- 
tributions, Nh = 1 and the proportionality factor is —28/79. With the entropy density 
s = 7.04n 7 in terms of the photon density, the baryon asymmetry of the Universe, defined 
by the ratio r\ B of the net baryon number to the photon number, is given in terms of the 
lepton asymmetry and the washout parameter Ki by 

v sm = ^£ _ _ Q216ki€i (29) 

Successful leptogenesis will require that the final result for i]b should lie in the observed 
range [23] 

T] B = (6.15 ±0.25) x 1(T 10 . (30) 

In the case where supersymmetric particles are considered to contribute to the decays, 
entropy and sphaleron interactions, the following modifications are in order: (a) The presence 
of sleptons and higgsinos in the self-energy and vertex loops will double the interference terms 
without affecting the tree-level decay rates, so the asymmetry is doubled, (b) The presence 
of sleptons and higgsinos in the decay products of the right-handed neutrinos will double 
the decay rates without affecting the asymmetry, (c) With the sneutrino counterparts of the 
decaying heavy neutrinos taken into account, the lepton asymmetry is further doubled, (d) 
The value of the effective number of relativistic degrees of freedom is now g* = 228.75. (e) 
In Eq. (28) we must take N H = 2 for the two Higgs doublets present in the supersymmetric 
case. 

The result is that the lepton asymmetry is replaced by a factor of 4 times the value 
of ei given in Eq. (26), but since the widths of the decaying states are doubled, the mass 
separation should be taken twice as large in the computation of the asymmetry. The washout 
factor as determined by the Boltzmann equations is also affected, but in a very complicated 
fashion with the supersymmetric particles present [24]. We shall assume that the expression 
in Eq. (24) is still a good approximation although this is open to question. The appropriate 
expressions now read 

SUSY 

~ 1.62 x 10 -2 Kiei, (31) 

SUSY 

~ -5.62 x 10 _3 Kiei, (32) 
v susy _ _o.0396«iei. (33) 

We now apply the above framework to the 50(10) model of [9,16] to determine just 
how successful leptogenesis would be in inducing baryogenesis given the mass matrices and 
parameters that lead to the observed neutrino masses and mixings. In this model, most of 
the complex phases can be rotated away from the Dirac mass matrices M Uy M D , M L , and 
Mtv; in fact all but two can, which were called a and <fi in [9]. The phase a gets set to zero 
by fitting the quark masses which is why it does not show up in Eq. (2). The phase is 
more important and is responsible for the CP violating phase Sqkm i n the CKM matrix. 
In order to get substantial leptogenesis in this model there must be a large phase in M R . 
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Now, in [16] Mr was taken to be real, as it was in in Sect. II, simply for convenience of 
analysis, as the interest there was only atmospheric and solar neutrino oscillations rather 
than leptonic CP-violating effects. An example of the fit with real Mr is given in Table 1 as 
case (1-0). Note that with the type I seesaw mechanism only four parameters, a, b, c and 
Ar, are required in addition to the eight present in the four Dirac mass matrices to obtain 
excellent agreement with the neutrino mass and mixing data. However, there is no reason 
for Mr to be real, and so here we shall investigate complex values. 
A simple choice of complex parameters for Mr is 

a = 1.2 -0.45i, b = 2.0, , . 

c = 2.0, A R = 2.65 x 10 14 GeV. [ ' 

This we call case (1-1), and the resulting observables are shown also in Table 1. Note 
that while the neutrino mixing parameters are all within the present acceptable range, the 
baryon asymmetry t]b is much too small, as it is almost four orders of magnitude below 
the observed value. This can be traced to the fact that the relative separation of the two 
quasi-degenerate right-handed neutrino masses is too large compared with the widths of 
their levels, (M 2 - M 1 )/T 1 = 34.2. 

The best results we have been able to obtain in the model arise if we allow small non-zero 
values in the 12, 13, 21, and 31 entries of the Dirac neutrino matrix. Consider, for example, 
the choice 

a = 0.25 + 0.15i, 6=1.2 + 0.9i, 

c = 0.25 + 0.25i, A R = 2.9 x 10 14 GeV, (35) 

(M N ) 12 = (M N ) 21 = -0.65 x lO" 5 , (M N ) 13 = (M N ) 31 = -1.0 x lO" 5 , 

which we call case (1-2). We have checked that the introduction of these non-zero values into 
the Dirac neutrino mass matrix, and likewise for the up quark mass matrix, does not destroy 
the good agreement for the quark masses and CKM mixings. We see from Table 1 that the 
neutrino mixing parameters for case (1-2) are also in very good agreement with the present 
known data. The splitting of two quasi-degenerate right-handed neutrino masses has now 
been reduced to half the width of either state, which maximizes the resonance enhancement. 
Note that the relative CP-parity of the two lightest neutrino states is no longer opposite 
but differs from that by about 4°, a necessary condition to get satisfactory leptogenesis. 
Moreover, for this type I seesaw model, the Dirac phase 5c p is large and much closer to 
maximal than in case (1-1). But the baryon asymmetry is r] B = 2.2 x 10 -10 , which falls short 
of the observed value by a factor of three. The biggest improvement over case (1-1) occurs in 
the lepton asymmetry which has improved by three orders of magnitude, while the washout 
factor is only slightly larger. 

The numbers given for the three cases discussed so far for cases (1-0), (1-1), and (1-2) all 
left out the contributions from the supersymmetric particles. In case (1-3) of Table 1, we 
have taken those contributions into account. For this case we take 

a = 0.2 + 0.H, 6 = 1.15 + 0.9i, 

c = 0.25 + 0.3i, A R = 2.84 x 10 14 GeV, (36) 

{M N ) 12 = (M N ) 21 = -0.65 x lO" 5 , (M N ) 13 = (M N ) 31 = -1.05 x 10" 5 , 
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We see from Table 1 that although tjb given by Eq. (33) in the SUSY case appears to be 
twice as large as in the SM calculated from Eq. (19) e\ is only half as large. The net 

baryon asymmetries are thus nearly equal for cases (1-2) and (1-3), provided one is justified in 
using the same washout formula in both cases. It thus appears that, in this realistic model 
of quark and lepton masses and mixings, obtaining sufficient baryon asymmetry through 
thermal leptogenesis is somewhat problematic. 



IV. LEPTOGENESIS WITH TYPE III SEESAW 



In a recent paper [12] a new type of seesaw mechanism was proposed for light neutrino 
masses that can be implemented in grand unified theories based on SO (10) or larger groups. 
This was called the type III see-saw mechanism and can be implemented in the model of 
[9,16], which we have been considering. Indeed, in [13] it was argued that when implemented 
in this model, the type III see-saw can give both realistic neutrino masses and mixings and 
sufficiently large leptogenesis without fine-tuning of the forms of the matrices. Here we will 
look at type III leptogenesis in more detail, giving numerical examples. 

The type III seesaw involves introducing, in addition to the left- and right-handed neu- 
trinos (z/j, iVj) contained in the 16j, three SO (lO)-singlet neutrinos 1$. Thus the neutrino 
mass matrix is not 6x6 but 9x9, and is given by 



(37) 



where the ^ C 5(16$) are the usual left-handed neutrinos, the iV? = 1(16$) are conjugates 
of the usual right-handed neutrinos, and the S$ = 1$ are the SO (lO)-singlet neutrinos. The 
index % runs over families. The submatrix represents the usual Dirac mass matrix 
contribution involving the doublet and singlet neutrinos in the 16j's as in the previous Sect. 
The new terms in the third row and third column arise from couplings involving the three 
singlet neutrinos Si as given by the additional contributions to the Yukawa superpotential: 





( (M N ) i:j 


F> jU 


\ 


( V 3 


\ 


W neut = (vi,N?,Si) 


(M&ij 


FijQ 




N j 






{ (F' T ) l]U (F^Q 


Mij 


J 




J 



w. 



RHv 



Mijhlj. 



(38) 



The superscript a distinguishes the different Higgs 16# representations, if there are more 
than one. The F' matrix in Eq. (37) arises when at least one of the Higgs fields, 5(16#), 
gets an electroweak VEV, u a , in the 5 SU(5) subgroup direction. Likewise, the F matrix 
appears when at least one of the 1(16#) fields get a superheavy VEV, Vl a) in the SU(5) 
singlet direction. More explicitly, we mean 



T F a Q 



FijQ 



(39) 



With the submatrix Ai also superheavy, W neu t can be diagonalized to yield the 3x3 matrix 
in the light left-handed neutrino vv block 



M„ 



u 



-M.xM^M.i- - (M N H + H 1 M^) — , 



(40) 
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where 



M R = (FQ)M-\F T Q), 
H = (F'F-y. 



(41) 



The first term in Eq. (40) is the usual type I seesaw contribution. The second term is 
the new type III seesaw contribution and arises when the 16^ fields develop an electroweak 
VEV, so F' 7^ 0. If the elements of the matrix M. are small compared to those of FVt, then 
it is easy to see from Eqs. (40) and (41) that the type I contribution becomes negligible 
compared to the type III contribution. In the limit that M. = one sees from Eq. (37) 
that the superheavy neutrinos have simply the mass term FijQ(NfSj). That is, the Nf and 
Si pair up to form three Dirac neutrinos. On the other hand, if Ai is small (compared to 
FQ) but not zero, then these three Dirac neutrinos get slightly split into six eigenstates 
forming three nearly degenerate pseudo-Dirac neutrinos. It is this fact that can be exploited 
to enhance leptogenesis. 

In Ref. [13] it was assumed that the matrices F, F' and Ai in the original flavor basis 
all have elements of the order 



(42) 



where A = rj/e = 4.1 x 10~ 5 . This form is suggested by that of M/v in Eq. (2), where the 
11 element is much smaller than the other non-zero elements, possibly due to an Abelian 
flavor symmetry. It is convenient to go to a basis where F is diagonalized. This is done by 
a biunitary transformation. We indicate quantities in this basis by the tilde symbol. These 
are related to those in the original flavor basis by 




F = U T FV, N c = UN C , S = VS. 



Then 



FQ = 





V o 





F 2 




where the Fj are of order unity. And 

M = | Xgi2 922 923 
Xgi3 g23 <?33 



\ 




(M x 





\ 







n = 





M 2 







F 3 J 









mJ 












Xf[ 2 




5, 


F'u 


■( 


A/21 


f22 


/23 








A/31 


/32 


/33 



(43) 



(44) 



(45) 



Again we assume //•, ~ 1. Here M s <C fl in order to obtain three superheavy quasi-Dirac 
neutrino pairs. 

Finally we transform the Dirac neutrino mass matrix according to M N = M N U. Because 
of the assumed form of F, the matrix U has the form 



13 



u 



I Mil Xu 12 XU13 \ 

Xu 2 i u 2 2 u 2 3 
V Aw 3 i U 32 M33 / 



where ~ 1 and by unitarity u n ~ I. The transformed Dirac neutrino matrix 
becomes 



M 



N 



1 rju n ^A-ui2 r]Xu 13 
eAw 3 i eu 32 eu 33 | m v = Ymu. 

\ Xu 31 U 32 U 33 



The type III seesaw mechanism given in Eq. (40) with dominance of the second term 
yields for the light left-handed neutrino mass matrix (recall that r]/X = e) 



(46) 
then 

(47) 
then 



M v = — 



'111/2 



2e Ej 

M 3j/ 2j 



/ mi/ 21 



"3j/ 3j -"2j/ 2 ,, 



+ e 



«3j/3 7 --«2j7 2 



2i 




ft 



(48) 



This clearly has a well-defined hierarchical form which is similar to the corresponding M v 
determined in Sect. III. for the type I seesaw mechanism. 

Leptogenesis is almost exclusively produced by the decays of the lightest pair of the six 
superheavy neutrinos. Neglecting, as is justified, the mixing of the lightest pair of superheavy 
neutrinos with the two heavier pairs of superheavy neutrinos, we find the effective two-by-two 
mass matrix for the lightest pair to be 



(JV?,5i) A 2 



Fifi 
Fifi gil M s 



Si 



(49) 



If, as we assume, Ms <C ft, these form an almost degenerate pseudo-Dirac pair, or equiv- 
alently two Majorana neutrinos with nearly equal and opposite masses. These Majorana 
neutrinos are JVi± = [N{±Si)/y/2, with masses M 1± = ±M l + \Mn = X 2 (±F 1 Vt+\g ll M s ). 
These can decay into light neutrino plus Higgs boson via the term Y i ±(Ni±h / i )H, where 



y 1± = (y i1 ±f; 1 )/V2t 



Mu (Y =F' 



2. 



(50) 



Here Y is the Dirac Yukawa coupling matrix given in Eq. (47). 

It is straightforward to show that the lepton asymmetry per decay produced by the 
decays of iV 1± is given by [25,5] 



1 imv ; iv; 

47rE,-[|^ + l 2 + l^- 



-J{Ml + /Ml 



(51) 



where f(Mf + /Mf_) comes from the absorptive part of the decay amplitude of N± and 
was given earlier in Eq. (17). For the application here f(Mf + /Mf_) = —Mi/2M n = 
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[Fi/2gi)(Vt/ M s ) which can be large if M s <C 0,. The expression for /(Mf + /Mf_) given 



\M_ 



Jvi\i is larger than half the widths 

2 



above is only valid when the mass splitting \M + 

of the N 1± , which from Eq. (20) are given by r± = £ fc |Yfc±| 2 . From Eqs. (45), (47), 
and (50), one sees that T\± ~ \ 2 Mi/8ti. As we shall see with our numerical examples, this 
condition that the splitting of Ni± be comparable to or greater than their widths, is easily 
satisfied. 

Making use of Eqs. (50) and (51) one obtains 



1 EjdYj! 


2 _ 




IXE t y;A) 




2 + 


3i 


2 ) 



f(Ml/Ml), 



(52) 



This can be evaluated in terms of the parameters of the model using Eqs. (45) and (47), 
giving 



A 2 

47T 



(l«3l| 


2 - l/£il 2 )M«5i&) 




«31 


I 2 + \& 


I 2 + \fL 


2 



/(M 2 + /M 2 _). 



(53) 



The washout parameter is given approximately as before by Eq. (24) where now 

8™ 2 IV ± 



mi 



M 2 m± 



A 2 ^(Kil 2 + l^l 2 + l^l 2 ). 



(54) 



The lepton asymmetry is then translated into the baryon asymmetry for the Universe by 
the same formulas which appeared earlier in Sect. III. 

We now wish to give some numerical examples for this type III seesaw mechanism. The 
first requirement, of course, is that we choose values of the parameters that reproduce the 
neutrino mass and mixing data. One could search over the whole space of parameters, 
but this is a cumbersome task, as in the type III mechanism there are many parameters 
involved in the sector of superheavy singlet neutrinos, including the parameters and f[y 
A more convenient approach is to find values of u.ij and f-j that make the matrix M v (given 
approximately in Eq. (48)) have a form close to that shown in Eq. (6), since we already 
know that that form can reproduce the neutrino mass and mixing data for suitable a, b, and 
c. First let us choose values of that are simple and such that U in Eq. (46) is unitary. 
They can not be too simple, i.e., have too many zeros, or else there will not be sufficient 
leptogenesis. We choose the following form for simplicity: 



U 



( 1 

-A(l + V2)i 



-A(l + y/2)i 
1/V2 
i/V2 



A \ 
1/V2J 



(55) 



where A = rj/e = 4.1 x 10 5 as before. We can make M v very close numerically to the form 
in Eq. (6) by taking 
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fll 

/21 

^23 
/3I 

^33 



0, with f[ 2 and f[ 3 arbitrary, 
-FjQ-a), 



V2 Fi 



F 2 /(b-a), 



1 " 272* (W + i73 (t?) + ( & " ] - «), 



(56) 



0, 
1 

2^/2 
1 

2\/2 



-i + 2(2 + >/2) 
1 + 2y/2&] F 3 , 



:1m. 



For simplicity we take the M. matrix of Eq. (45) to be diagonal and set gu = g 2 2 = #33 = 1- 
Three examples are illustrated in Table 2, with the values of f 21 , f 22 , f'23 determined from Eq. 
(56) by setting a, b, c set equal to 1.2 - 0.45i, 2.0, 2.0, respectively, as in case (1-1) of Table 
1; in addition Q appearing in the type III seesaw is set equal to A# = 2.7 x 10 14 appearing 
in the case (1-1). In these examples, F 1: F 2) F 3 and M s are allowed to vary. We have 
considered only the contributions to leptogenesis from the SM particles, since the previous 
results indicated little difference with the inclusion of the supersymmetric contributions, 
given the uncertainty in the corresponding washout factor. 

By using the type III formulas beginning with the transformed version of the light neu- 
trino mass matrix given in Eq. (40), we clearly reproduce the neutrino mass and mixing 
results obtained with the type I seesaw mechanism in Table 1 for Case (1-1), where no mod- 
ification of the original Dirac neutrino mass matrix was introduced. In fact, these results 
are independent of the actual values taken for Fi, F 2 , F 3 and Ms, and their values are not 
repeated in Table 2. In particular, the Dirac CP- violating phase remains far from maximal 
at 172°. On the other hand, one sees that successful leptogenesis and observed baryogen- 
esis are easily obtained in all three examples for those four input parameters chosen. The 
masses of the lightest pseudo-Dirac pair of heavy neutrinos range from 4.5 x 10 4 GeV up to 
4.5 x 10 6 GeV. For this mass range, overproduction of gravitinos is not an issue [26]. The 
mass separation of the two Majorana neutrinos in the lightest pair is well above their decay 
widths. Hence we see that the type III seesaw mechanism has a much easier time simultane- 
ously giving realistic light neutrino masses and mixings and successful leptogenesis. In the 
realistic model we have examined, we see that it fails by at least a factor of three to give 
enough leptogenesis with the conventional seesaw mechanism. By contrast, with the type 
III seesaw mechanism as we have seen, the light neutrino mixing issue is decoupled from 
the issue of leptogenesis; moreover, there are ready-made pairs of nearly degenerate heavy 
neutrinos. Consequently, there is no difficulty in the same realistic model getting successful 
leptogenesis with the type III seesaw mechanism. 



V. CONCLUSIONS 

In this paper we have explored the issue of resonant leptogenesis in a very predictive 
5*0(10) grand unified model which leads to reliable numerical results. Phenomenological 
studies [6] by Ellis, Raidal, and Yanagida, as well as by Akhmedov, Frigerio, and Smirnov 
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have previously suggested that successful leptogenesis can be easily obtained with resonant 
enhancement, though they had not looked at any specific realistic models in detail to be 
sure that both the leptogenesis and neutrino mass and mixing results can be simultaneously 
satisfied. In order to obtain the desired LMA solar neutrino solution in the model con- 
sidered here, it was found that near degeneracy of the lightest two right-handed neutrinos 
was required. This feature neatly favors the resonant enhancement scenario for leptogenesis 
which can survive washout and be converted into the observed baryon excess by sphaleron 
interactions in thermal equilibrium above the critical electroweak symmetry-breaking tem- 
perature. 

We have first studied this problem in the conventional type I seesaw framework. The 
seven model parameters, including one complex one, in the Dirac matrices are fixed by 
the charged lepton and quark mass and mixing data. The lopsided nature of the down 
quark and charged lepton mass matrices neatly accounts for both the small value of the 
V c b quark mixing parameter and the near maximal mixing of the mu- and tau-neutrinos. 
As originally proposed, the right-handed Majorana neutrino mass matrix depended on just 
four real parameters which, in conjunction with the Dirac neutrino and charged lepton mass 
matrices, then leads to a sub-maximal mixing of the solar neutrinos but no leptogenesis. 

By allowing three of the four M R parameters to be complex and introducing two addi- 
tional very small parameters into the Dirac neutrino mass matrix, we were able to achieve 
a sizeable amount of leptogenesis; however, it falls short of the observed value by a factor of 
three. The limiting factor is the requirement that the heavy neutrino M\ — M2 mass separa- 
tion must be comparable to or larger than half the decay width of either neutrino. For this 
application the two masses are found to be of the order of 3 x 10 8 GeV and separated by 600 
GeV. In the supergravity scenario of SUSY breaking, however, the gravitino problem which 
requires an upper bound on the reheating temperature of T R < 10 7 GeV at best makes this 
solution somewhat problematic. This can be alleviated, however, if SUSY breaking occurs 
via the gauge-mediated scenario. 

We then considered a type III seesaw mechanism by which three 5*0(10) singlet neutri- 
nos are added to the spectrum, so the neutrino matrix is 9 x 9. With the singlet neutrinos 
at some large intermediate scale, the effective "double seesaw" results in three pairs of 
quasi-Dirac neutrinos in place of three right-handed Majorana neutrinos. In this scenario 
resonant-enhanced leptogenesis is achieved by the lightest pair of quasi-Dirac neutrinos, for 
which their mass separation is typically one hundred times the widths of the two states. 
Hence successful leptogenesis is easily obtained. With this type III seesaw the leptogenesis 
is completely decoupled from the neutrino mass and mixing issues, so the good agreement of 
the latter with present-day observations is preserved as observed earlier by us in [13]. Here 
we have seen that with one choice of parameters the masses of the lightest neutrino pair can 
be of the order of 5 x 10 4 GeV, so the gravitino problem is no longer an issue. The only 
drawback for this type III seesaw mechanism is the appearance of a large number of new 
parameters needed to specify the expanded 9x9 neutrino matrix. 

The research of SMB was supported in part by the Department of Energy under contract 
No. DE-FG02-91ER-40626. One of us (CHA) thanks the Theory Group at Fermilab for 
its kind hospitality. Fermilab is operated by Universities Research Association Inc. under 
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TABLES 

TABLE I. Type I seesaw results for the four cases considered in the text. 





Case (1-0): (SM) 


Case (LI): (SM) 


Case (1-2): (SM) 


Case (1-3): (SUSY) 


Input: 


a 


1.0 


1.20 - 0.45i 


0.25 + 0.15i 


0.20 + O.lOi 




b 


2.0 


2.0 


1.20 + 0.90i 


1.15 + 0.90i 




c 


2.0 


2.0 


0.25 + 0.25i 


0.25 + 0.3i 




A R (GeV) 


2.65 x 10 14 


2.70 x 10 14 


2.90 x 10 14 


2.84 x 10 14 




( 7t/f ^ ( 71 xT "\ 
(M N )i2, {M N )21 


0.0 


0.0 


-0.65 x 10" 5 


-0.65 x 10~ 5 




\M N Ji3, (M N )3i 


0.0 


0.0 


-1.0 x 10~ 5 


-1.05 x 10 -5 


Uutput. 


Nl\ ) 


2.31 x 10 8 


2.16 x 10 8 


3.06 x 10 8 


3.08 x 10 8 




1V12 [KjQ V ) 


2.31 x 10 8 


2.16 x 10 8 


3.06 x 10 8 


3.08 x 10 8 




A/f„ iTdV 1 ! 
(KjrQ V ) 


2.71 x 10 14 


2.76 x 10 14 


2.90 x 10 14 


2.86 x 10 14 






1.21 x 10~ 4 


1.31 x 10~ 4 


1.96 x 10~ 6 


3.87 x 10~ 6 




1 l/Mi 


3.83 x 10 b 


3.83 x 10~ 6 


3.84 x 10 b 


3.83 x 10~ 6 




vtl\ (meV) 


5.1 


5.3 


2.8 


2.8 




TYI2 (meV) 


9.1 


9.8 


8.8 


8.8 




7713 (meV) 


52.0 


52.0 


51.0 


52.0 




Am 32 (ev ) 


2.6 x 10~ 3 


2.6 x 10~ 3 


2.5 x 10~ 3 


2.6 x 10~ 3 




Am 21 (ev ) 


5.6 x 10~ 5 


6.9 x 10~ 5 


6.9 x 10~ 5 


6.9 x 10~ 5 




Sin ^Patm. 


0.991 


0.988 


0.979 


0.981 




„; T1 2 0/1 

sm zy so i 


0.87 


0.86 


0.80 


0.80 




xan (7^2 


0.48 


0.46 


0.39 


0.39 




Sin Zfi3 


0.0012 


0.0006 


0.0012 


0.0014 




IT/ „| 


0.017 


0.012 


0.018 


0.019 






177° 


172° 


76° 


71° 




Xl, X2 


-180°, 90° 


24°, 117° 


102°, -170° 


102°, -170° 




(m ee ) (meV) 


0.57 


0.58 


0.44 


0.44 




ei 


0.0 


-6.4 x 10~ 7 


-1.2 x 10" 3 


-6.1 x 10~ 4 




mi (eV) 


12.6 


13.5 


9.5 


9.5 






6.2 x 10" 6 


5.8 x 10~ 6 


8.3 x 10~ 6 


8.4 x 10~ 6 






0.0 


0.80 x 10~ 13 


2.2 x 10~ 10 


2.0 x 10~ 10 
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TABLE II. Type III seesaw results for three examples patterned after case (1-1) in Table 1. 



The light neutrino mass and mixings results obtained are identical to those 


in case (1-1) of Table 1 


and are not repeated here. 










Case (III-l): (SM) 


Case (III-2): (SM) 


Case (III-3): (SM) 


Input: Q (GeV) 


2.7 x 10 14 


2.7 x 10 14 


2.7 x 10 14 


Fi 


1.0 


10. 


0.1 


F 2 


1.0 


0.1 


0.1 


F 3 


1.0 


1.0 


1.0 


M s 


4.3 x 10 5 


8.5 x 10 8 


1.0 x 10 5 


fu/Fi 


0.0 


0.0 


0.0 


/12/-P2 


1.0 


1.0 


1.0 


/13/-P3 


1.0 


1.0 


1.0 


/21/F1 


-0.950 + 0.534i 


-9.496 + 5.3411 


-0.095 + 0.053i 


/22/-F2 


-2.279 - 1.537i 


-0.227 - 0.1541 


-0.228 - 0,154i 


/23/-P3 


-0.194 + 1.523i 


-0.194 + 1.523i 


-0.194 + 1.523i 


fsi/F! 


0.0 


0.0 


0.0 


/32/-F2 


-0.354i 


-0.035i 


-0.035i 


/33/ -^3 


0.354 


0.354 


0.354 


511 


1.0 


1.0 


1.0 


522 


1.0 


1.0 


1.0 


533 


1.0 


1.0 


1.0 
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FIGURES 





Figure 1: Tree-level and one-loop Feynman diagrams in the computation of 
the lepton asymmetry. A similar set of diagrams exists for the decays into 
the antineutrino channels. 
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